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Abstract 

We investigate the ferromagnetic-glassy transitions which separate the low-temperature ferro- 
magnetic and spin-glass phases in the temperature-disorder phase diagram of three-dimensional 
Ising spin-glass models. For this purpose, we consider the cubic-lattice ±J (Edwards-Anderson) 
Ising model with bond distribution P(J) = pS(J — 1) + (1 — p)S(J + 1), and present a numerical 
Monte Carlo study of the critical behavior along the line that marks the onset of ferromagnetism. 

o' 

O . The finite-size scaling analysis of the Monte Carlo data shows that the ferromagnetic-glassy 

transition line is slightly reentrant. As a consequence, for an interval of the disorder parameter p, 

> 

l/~) ' around p ~ 0.77, the system presents a low-temperature glassy phase, an intermediate ferromag- 

o" 

netic phase, and a high-temperature paramagnetic phase. Along the ferromagnetic-glassy transition 

en 

i^ , line magnetic correlations show a universal critical behavior with critical exponents v = 0.96(2) 



CO 



o 



and n = —0.39(2). The hyperscaling relation fi/v = (1 + n)/2 is satisfied at the transitions, so that 
f3/v = 0.305(10). This magnetic critical behavior represents a new universality class for ferromag- 



S^ netic transitions in Ising-like disordered systems. Overlap correlations are apparently not critical 

h ; 

and show a smooth behavior across the transition. 
PACS numbers: 75.50.Lk,05.70.Fh,64.60.F-,05.10.Ln 



I. INTRODUCTION 

Spin glass models are simplified, although still quite complex, models retaining the main 
features of physical systems which show glassy behavior in some region of their phase di- 
agram. They may be considered as theoretical laboratories where the combined effects of 
disorder and frustration can be investigated. Their phase diagram and critical behavior 
can be used to interpret the experimental results for complex materials. Ising-like spin 
glasses, such as the ±J Ising model, 1 model disordered uniaxial magnetic materials char- 
acterized by random ferromagnetic and antiferromagnetic short-ranged interactions, such 
as Fei-^Mn^TiOs and Eui-^Ba^MnOs; see, e.g., Refs. |2|-|4|. The random nature of the 
short-ranged interactions is mimicked by nearest-neighbor random bonds. 

Three-dimensional (3D) Ising spin glasses have been widely investigated. At low temper- 
atures they present ferromagnetic and glassy phases, depending on the amount of frustra- 
tion. The critical behaviors along the finite-temperature paramagnetic-ferromagnetic and 
paramagnetic-glassy (PG) transition lines have been accurately studied.- - — On the other 
hand, the low-temperature behavior, in particular the nature of the glassy phase and of the 
boundary between the ferromagnetic and glassy phases, is still debated. 

In this paper we focus on the low-temperature transition line which separates the ferro- 
magnetic phase, characterized by a nonzero magnetization, and the spin-glass (glassy) phase 
in which the magnetization vanishes while the overlap expectation value remains nonzero. 
We consider the 3D ± J Ising model, defined by the Hamiltoniani 

H — — y j Jxy&x&y, (1) 

(xy) 

where a x = ±1, the sum is over the nearest-neighbor sites of a cubic lattice, and the exchange 
interactions J xy are uncorrelated quenched random variables with probability distribution 

P{J xy ) = p5(J xy - 1) + (1 - p)S(J xy + 1). (2) 

The usual bimodal Ising spin glass model, for which [J xy \ = (brackets indicate the average 
over the disorder distribution), corresponds top = 1/2. For p ^ 1/2 we have [J xy ] = 2p—l ^ 
0, and ferromagnetic (or antiferromagnetic) configurations are energetically favored. 

The phase diagram of the cubic-lattice ±J Ising model is sketched in Fig. [TJ We only 
consider p > 1/2 because of the symmetry p — > 1 — p. While the high-temperature phase is 




FIG. 1: (Color online) Temperature-disorder phase diagram of the 3D ±J Ising model. The phase 
diagram is symmetric under p — ^ 1 — p (but for small values of p the system is antiferromagnetic) . 



always paramagnetic (P), at low temperatures there is a ferromagnetic (F) phase for small 
frustration, i.e., small values of 1 — p, and a glassy (G) phase with vanishing magnetization 
for sufficiently large frustration. In Fig. [1] we do not report any low-temperature mixed phase 
with simultaneous glassy and ferromagnetic behavior as found in mean-field models^, for 
which, at present, there is no evidence.— 1 ^ The different phases are separated by transition 
lines belonging to different universality classes. They meet at a magnetic-glassy multicritical 
point M located along the so-called Nishimori line^ 1 ^- 2/T = ln[p/(l — p)], where the 
magnetic and the overlap two-point correlation functions are equal. Scaling arguments^ 1 ^ 
show that the transition lines must be all parallel to the T axis at the multicritical point M. 
The paramagnetic-ferromagnetic (PF) transition line starts at the Ising transition of 
the pure system at p — 1, at 30 Tj s = 4.5115232(16), with a correlation-length exponent 
z/fe « 0.6301 (u = 0.63012(16) from Ref.blland v = 0.63002(10) from Ref.boJ). Along the PF 
line the magnetic critical behavior is universal,- and belongs to the randomly-dilute Ising uni- 
versality class,— >22 characterized by the correlation-length critical exponent z/pp = 0.683(2). 
It extends up to the multicritical point M, located alr^ pm = 0.76820(4), T M = 1.6692(3), 
whose multicritical behavior is characterized by two even relevant renormalizat ion-group 
(RG) perturbations with RG dimensions y\ = 1.02(5) and y^ = 0.61(2). The paramagnetic- 
glassy (PG) transition line runs from M to the finite-temperature transition at p — 1/2, 
at 6 Tb = 1.11(1). The glassy critical behavior is universal along the PG line;- the overlap 



correlation-length exponent is quite large,- - — upq = 2.45(15). Finally, (at least) another 
transition line is expected to separate the ferromagnetic and glassy phases. This is the 
ferromagnetic-glassy (FG) transition line that marks the onset of ferromagnetism and which 
runs from M down to the point D at T = 0. The nature and the general features of this tran- 
sition line in Ising spin glasses are not known. Beside a few numerical works at T = 0,— ^ 
this issue has never been investigated at finite temperature. 

An interesting issue concerning the FG transition line is whether it is reentrant, which 
would imply the existence of a range of values of p for which the glassy phase is separated 



from the paramagnetic phase by an intermediate ferromagnetic phase. As proved in Refs. 



id, 



171 . ferromagnetism can only exist in the region p > pm, which implies that p^ > pu- We 
also mention that, using entropic arguments applied to frustration, the FG phase boundary 
was argued to run parallel to the T axis,— ^ i.e., pn = pm for any T < Tm, with the 
critical behavior controlled by a T = percolation fixed point. 13 The FG transition was 



numerically investigated at T = in Ref. 



14 . obtaining the estimate j?d = 0.778(5) for the 



critical disorder, which is slightly larger than pu = 0.76820(4). Thus, it suggests a slightly 
reentrant FG transition line, although its apparent precision is not sufficient to exclude 

Pd = Pm- 

In this paper we study the nature of the FG transition. In particular, we investigate 
whether the magnetic variables show a continuous and universal critical behavior from M to 
D, and whether hyperscaling is violated as it occurs in some systems whose critical behavior 
is controlled by a zero-temperature fixed point, like the 3D random-field Ising model.— 

Note that we focus on the low-temperature ferromagnetic transition line, which marks 
the onset of ferromagnetism moving from the glassy phase with zero magnetization. There 
is also the possibility that a second low-temperature transition line exists for larger values 
of p. In this case there would be a mixed low-temperature phase, in which ferromagnetism 
and glassy order coexist. This occurs in mean-field models^ such as the infinite-range 
Sherrington-Kirkpatrick model.— However, numerical T = ground-state calculations in 
the 3D ± J Ising model on a cubic lattice— and in related models^ do not seem to show 
evidence of a mixed phase and are consistent with a unique transition. 

In this paper we present a Monte Carlo (MC) study of the critical behavior along the 
FG transition line. We perform simulations of finite systems defined on cubic lattices of size 
L < 20. A finite-size scaling (FSS) analysis of numerical data at T = 0.5 and T = 1 as 



a function of p shows that magnetic correlations undergo a continuous transition along the 
FG line. The critical behavior is universal, i.e., independent of T along the line. For the 
magnetic critical exponents we obtain v = 0.96(2) and 77 = —0.39(2). Moreover, hyperscaling 
is verified. The FG transition line turns out to be slightly reentrant. Indeed, we find 
p c = 0.7729(2) at T = 0.5 and p c = 0.7705(2) at T = 1, which are definitely larger than 
the disorder parameter pu = 0.76820(4) at the multicritical point. Therefore, for a small 
interval of the disorder parameter, around p ~ 0.77, the phase diagram presents three 
different phases: a low-temperature glassy phase, an intermediate ferromagnetic phase, and 
a high-temperature paramagnetic phase. 

Note that the critical behavior of the magnetic correlations along the FG transition line 
shows a new universality class of ferromagnetic transitions in Ising-like disordered systems, 
which differs from the randomly-dilute Ising universality class describing the critical behavior 
along the PF transition line, and from the random-field Ising universality class characterized 
by hyperscaling violation. 

The general features of the phase diagram presented in Fig. [I] should also characterize the 
temperature-disorder phase diagram of other 3D Ising spin glass models with tunable disor- 
der parameters. For example, one may consider models with Gaussian bond distributions, 

such as 

{Jxy ~ Jo) 



P(Jxy) ~ exp 



2a 



(3) 



where the parameters Jo and a control the amount of disorder (the pure ferromagnetic model 
corresponds to Jo > and a = 0). This distribution is also characterized by the presence of a 
Nishimori line T = a j Jo, where the magnetic and the overlap two-point correlation functions 
are equal. We also mention that an analogous temperature-disorder phase diagram, with 
three transition lines meeting at a multicritical point like Fig. (TJ is also found in 3D XY 
gauge glass models.— A similar phase diagram is also expected for other continuous spin 
glasses, like XY and Heisenberg spin glasses with bond distributions fl2]) or (jHJ). 

The paper is organized as follows. In Sec. HTIwe describe the MC simulations, and provide 
the definitions of the quantities we consider. Sec. IIHI presents the FSS analysis of the MC 
data, reporting the main results of the paper. Finally, in Sec. IIVI we draw our conclusions. 
In the appendix we report some details of the FSS analyses. 



II. MONTE CARLO SIMULATIONS AND OBSERVABLES 

In order to study the FG transition line, which connects points M and D in Fig. [1] we 
perform MC simulations of the ± J Ising model on cubic lattices of size L with periodic 
boundary conditions. We use the Metropolis algorithm, the random-exchange method, and 
multispin coding. Implementation details can be found in Ref. |6j. In the random-exchange 
simulations we consider Nt systems at the same value of p and at different temperatures in 
the range T max > T- h > T min , with T max > 2 and T min = 0.5. The value T max is chosen so 
that the thermalization at T max is sufficiently fast — typically we take T max > T M « 1.67— 
while the intermediate values Tj are chosen such that the acceptance probability for the 
temperature exchange is at least 10%. We require one of the Tj to be along the Nishimori 
line.— The results for this temperature value can be compared with the known exact results 
and thus provide a check of the MC code and the thermalization. Finally, one of the 
temperatures always corresponds to T = 1. The parameter Nt increases with L and varies 
from Nt = 5 for L = 4 to Nt = 19 for L = 20. Thermalization is checked by verifying 
that disorder averages are stable when increasing the number of MC steps for each disorder 
realization. We average over a large number N s of disorder samples: N s m 2 x 10 6 samples 
for L = 4, 6, 8, N s « 3 x 10 5 for L = 10, N s « 10 5 for L = 12, N s « 5 x 10 4 for L = 16, and 
N s « 5 x 10 3 for L = 20. 

The simulations are quite costly, because of the very slow dynamics for low temperatures. 
This makes the computational effort increase with a large power of the lattice size. In our 
range of values of L, the number of iterations which must be discarded for thermalization 
apparently increases as L s for our largest lattices (with an increasing trend with increasing 
L). Hence, taking into account the volume factor, the CPU time for each disorder realization 
apparently increases as L u (but we should warn that its large-L asymptotic behavior may 
be even worse). In total, simulations took approximately 40 years of CPU time on a single 
core of a recent standard commercial processor. 

We consider the magnetization and the magnetic correlation function defined as 

™ = ^[<ll>l>], (4) 



V 1 



X 



G(x) = l(a or x )], 
where the angular and the square brackets indicate the thermal and the quenched average 



over disorder, respectively. We define the magnetic susceptibility and the second-moment 
correlation length, respectively as 

X 

1 G(0) - G(q) 



z 



4sin 2 (g min /2) G{q) 



where q = (q m in, 0, 0), g min = 2n/L, and G(q) is the Fourier transform of G(x). Moreover, 
we consider the cumulants 

* e & (6) 

^ 22 = ""EST"' 

where 

***= <£>«)*>■ (7) 

a; 

At the critical point R% = £/L, C/4, and ^22 (in the following we call them phenomenological 
couplings and denote them by R) are expected to approach universal values in the large-L 
limit (within cubic L 3 systems with periodic boundary conditions). In the ferromagnetic 
phase we have U4 — > 1, U22 -^ 0, and R% — > 00, while in the glassy phase we expect R^ — > 0. 

We also define analogous quantities using the overlap variables q x = a x a x , where a x 

(2) 
and a x are two independent replicas corresponding to the same couplings J xy . In particular, 

we consider £„ and JJ1 defined by replacing the magnetic variables with the overlap variables 

in Eqs. © and ©. 

III. FINITE-SIZE SCALING ANALYSIS 

In this section we present a finite-size scaling (FSS) analysis of the MC data close to the 
FG transition line. We consider two values of the temperature, T = 0.5 and T — 1, below 
the temperature Tm = 1.6692(3) of the multicritical point M, and perform a FSS analysis 
as a function of p. 
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FIG. 2: (Color online) Estimates of R^ at T = 0.5. The vertical lines show the location of the 
multicritical point M: 1 - p M = 0.23180(4). 



A. Phenomenological couplings and universality 

To begin with, we analyze the data at T = 0.5. In Fig. [2] we show the MC estimates 
of i?£ as a function of 1 — p. Analogous plots are obtained for U4 and 1/22- The data 
for different lattice sizes clearly show crossing points, providing evidence for a continuous 
transition. They cluster at values of p which are definitely larger than pm, ruling out a 
vertical transition line from M to the T = axis. 

In the critical limit, the phenomenological couplings R scale as 



R = f R [(p- Pc )L^] 



(8) 



where we have neglected analytic and nonanalytic scaling corrections. Equivalently, one can 
test FSS by considering two different couplings R\ and R2. In the FSS limit R\ = F^i^), 
where the function F^R-z) is universal, i.e., identical in any model that belongs to a given 
universality class. Clear evidence of FSS is provided in Fig. |3l where the phenomenological 
couplings C/4 and U22 are reported versus R^ = £/L, The data appear to rapidly approach 
a nontrivial limit with increasing the lattice size. Scaling corrections are only visible in the 
case of U22, but they decrease with increasing L. 

In order to determine the critical parameter p c and the exponent v, we fit C/4, U22, and 
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FIG. 3: (Color online) C/ 4 (bottom) and U 2 2 (top) vs i? ? at T = 0.5. 



i?£ = £/L to Eq. (JHJ). Details are reported in App. IA II We obtain 

Pc (T = 0.5) = 0.7729(2), v = 0.96(2), (9) 

Rl = 0.764(6), Ul = 1.331(5), U* 2 = 0.305(2), (10) 

where R* = /r(0) is the value of the phenomenological coupling R at the critical point. 
Scaling corrections turn out to be small. 

An analogous FSS analysis can be performed at T = 1, with the purpose of checking 
universality, i.e., of determining whether all transitions along the FG line belong to the same 
universality class. For this purpose, we use the fact that, given any pair of RG invariant 
quantities R± and R 2 , the FSS function Ri = Fi 2 (i?2) is universal. In Fig. 0] we plot [7 4 
and U22 versus R^ for both T = 0.5 and T — 1. The plot of U4 provides good evidence of 
universality: all data fall onto a single curve with remarkable precision. The results for U22 



show instead significant scatter, but they are also consistent with universality if one takes 
into account scaling corrections: indeed, as L increases the data for T = 1 approach the 
T = 0.5 results. 

For a more quantitative check, we must explicitly take into account scaling corrections at 
T = 1, since they are significantly larger than those observed at T = 0.5. For instance, fits of 
the phenomenological couplings at T = 1 to Eq. (jSj) show a somewhat large x 2 /DOF (DOF 
is the number of degrees of freedom of the fit). Moreover, the estimates show systematic 
trends as the lattices with smaller values of L are discarded in the fit, see App. IA II for 
details. To include scaling corrections, we fit the data to 

R = f R [(p - Pc)L l ' v ] + L-"g R [(p - p c )l}l"\. (11) 

The smallest y 2 /DOF is obtained for 0.8 < u < 0.9. Correspondingly v = 0.91(3), in 
substantial agreement with the estimate (EJ). Also the estimates of R^, U£, and U% 2 , see 
App. IA 1L are in agreement with the estimates ( TTOl) at T = 0.5. Therefore, all results 
strongly support the universality of the critical behavior along the FG line. It is difficult 
to estimate reliably the exponent u from the data. It we assume universality and fit the 
results at T = 1 fixing v = 0.96(2), we obtain u = 0.95(10). Note that the fits of the data 
at T = 0.5 give much larger values for u, i.e., u> > 2, see App. IA II This is probably due 
to the fact that corrections with w«l have very small amplitudes at T = 0.5, so that we 
are simply measuring an effective exponent that mimicks the behavior of several correction 
terms. 

The FSS fits also provide estimates of p c at T — 1. We obtain 

p c (T = 1) = 0.7705(2). (12) 

Note that p c {T = 1) > pu ~ 0.7682, conferming the reentrant nature of the FG transition 
line. 

B. Magnetic susceptibility 

As discussed at length in Ref. |6|, in the critical limit the magnetic susceptibility scales as 

X (p,L) = u h (p) 2 L 2 -V x [(p- Pc )L^}, (13) 
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FIG. 4: (Color online) U4 (bottom) and U22 (top) vs R^ at T = 1 and at T = 0.5 (only data with 
L> 12). 
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FIG. 5: (Color online) x = X u h £ versus £/L for T = 1 and T = 0.5. 



11 



where Uhip) is related to the magnetic scaling field and is an analytic function of p (and also 
of the temperature). Fits of x at T = 1 and T = 0.5 are good (x 2 /DOF of order 1) if we 
include all data such that L > 6, provided that Uh(p) is taken into account (see App. IA2I 
for details). We end up with the final estimate 

rj = -0.39(2). (14) 

Since £/L is a function of (p — p c )L l l v in the FSS limit, see Eq. (JSJ), we can rewrite Eq. ( TT31) 

as 

X (p,L)=u h (p) 2 e- ri F x (Z/L). (15) 

The function F x (x) is universal apart from a multiplicative constant, which takes into ac- 
count the freedom in the normalization of the function Uh(p). In Fig. [5] we show the quantity 
X = X u h £~ 2 ' 39 for T = 1 and T = 0.5. For each temperature the function Uh{p) is deter- 
mined by fitting the susceptibility data to Eq. ( ITBT) . fixing rj = —0.39. Moreover, the scaling 
fields are normalized so that x(T = 1,L = 16) ~ x(T = 0.5, L = 16) for £/L ~ 0.8. If we 
discard the data with L = 4 and 8 at T = 0.5, all points fall on top of each other, confirming 
universality. 

C. Evidence of hyperscaling 

Since the FG transition line extends up to T = 0, hence the critical behavior may be 

controlled by a zero-temperature fixed point, hyperscaling might be violated, as it happens 

in the 3D random-field Ising model.— In order to check whether hyperscaling holds along 

the FG line, we consider the magnetization, which is expected to behave as m ~ L~^l v at 

the critical point, and the magnetic susceptibility, which scales as x ~ L 2 ~ v . If hyperscaling 

holds, (3 and r/ are related by 

g = d - 2+ ", (16) 

v 2 K J 

(in the present case d — 3), which guarantees that xjvn 2 scales as L d . In order to verify 
whether Eq. f)16p holds, we consider H = %/(m 2 L 3 ) and assume that it behaves as 

H = -£-~tff H [(p-p c )LV% (17) 
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FIG. 6: (Color online) Scaling behavior of the distribution of the thermal averages of the magne- 
tization, at T = 0.5 and p = p c = 0.7729. We set P/v = 0.305. 

If hyperscaling holds, ( vanishes. A FSS analysis of the data at T = 0.5 and T = 1 gives 
the rather stringent bound (details in App. IA3j) 

|C| < 0.01, (18) 

which allows us to conclude, quite confidently, that hyperscaling holds. If this the case, 
using estimates (TT4I) and © of r\ and u, we obtain 

/3/u= (l + rj)/2 = 0.305(10), [5 = 0.29(1). (19) 

As a further check, we consider the sample distribution P(mt) of the thermal averages of 
the magnetization 



™* s 7<ll>*l> ( 20 ) 

X 

at the critical point p = p c = 0.7729, T = 0.5, which is expected to behave asymptotically 



as 



P(m t ) « l/ /v V(L^ /v m t ). (21) 

In Fig. O we plot V(L^^ u m t ) using /3/u = 0.305. The data clearly show the expected 
scaling behavior. In conclusion, the numerical results do not show evidence of hyperscaling 
violations in the critical behavior of magnetic correlations. 

Our data for H(p, L) can also be used to provide further evidence of universality. Indeed, 
if we use the fact that £/L is a function of (p — p c )L 1 / u , see Eq. (jSj), we can rewrite Eq. ([Tj 
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FIG. 7: (Color online) iJ = x/{m 2 L 3 ) versus £/L for T = 1 and T = 0.5. 

for ( = as 

H(p,L) = F H (£/L) + 0(L-"), (22) 

where Fh(x) should be the same at T = 0.5 and at T = 1 if all transitions along the FG 
transition line belong to the same universality class. The plot of the data, see Fig. [7J clearly 
confirms universality: all points fall onto a single curve. 

D. Overlap correlations 

In our numerical study we also consider quantities involving the overlap variables, such 
as £, /L and U%, defined at the end of Sec. [Til In Fig. [8] we show MC data up to L = 12 
(since their computation turned out to be significantly more demanding, we restricted the 
measurements for the lattices L = 16, 20 to the magnetic correlations). Unlike the magnetic 
quantities, the overlap data do not show crossings in the interval of p we have investigated. 
Apparently U% decreases continuously, while £ /L increases as L — > oo. This may reflect 
the fact that the FG transition line separates two ordered phases with respect to the overlap 
variables. Note that the differences between data at the same p and T and at different values 
of L decrease as 1 — p increases. Hence, if there is a line in the (p, T) plane where the overlap 
variables show crossings, it must be such that 1 — p > 0.234, i.e., it must lie in the region 
V < Pm, where no ferromagnetism is possible. 
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FIG. 8: (Color online) Estimates of £ /L (bottom) and U% (top), denned in terms of the overlap 
variables, at T = 0.5 and T = 1. 

IV. CONCLUSIONS 

We investigate the critical behavior along the ferromagnetic-glassy transition line of the 
T-p phase diagram of the cubic-lattice ±J (Edwards- Anderson) Ising model, cf. Eq. (PQ), 
which marks the low-temperature boundary between the ferromagnetic phase and the glassy 
phase where the magnetization vanishes, i.e., the transition line that runs from M down to 
the point D at T = in Fig. [Q 

We present a numerical study based on MC simulations of systems of size up to L = 20, 
obtaining MC estimates of several quantities at T = 0.5 and T = 1 (which are well below 
the temperature Tm = 1.6692(3) of the multicritical point M) as a function of the disorder 
parameter p. The results of the FSS analyses are consistent with the two continuous magnetic 
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transitions belonging to the same universality class. The corresponding critical exponents 
are v = 0.96(2) and 77 = —0.39(2). Since the critical line extends up to T = 0, the 
critical behavior may be controlled by a zero-temperature fixed point. Correspondingly, it 
is possible to have hyperscaling violations, as it occurs in the 3D random-field Ising model. 
Our MC results show that the hyperscaling relation ft/u = (1 + f])/2 is satisfied, so that 
fi/v = 0.305(10) and (3 = 0.29(1). The FSS results provide a robust evidence of a universal 
magnetic critical behavior along the FG transition line. A reasonable hypothesis is that also 
the zero-temperature transition belongs to the same universality class. This is supported 
by the available numerical data at T = 0. The numerical study of Ref. |l4J for the ±J 
Ising model at T = 0, using lattice sizes up to L = 14, provided evidence of a magnetic 
transition at pn = 0.778(5), with critical exponents v = 1.3(3) and /3 = 0.2(1). Numerical 
analyses^ for other Ising spin-glass models at T = give consistent values of the critical 
exponents, v = 0.9(2) and (3 = 0.3(1) using data up to L = 12. These estimates are 
substantially consistent with our results along the FG transition line, supporting a universal 
critical behavior along the FG transition from the multicritical point M down to the T = 
axis. 

We also investigate the behavior of overlap correlations. They do not appear to be critical 
and show an apparently smooth behavior across the FG transition. Our numerical results 
do not show evidence of other transitions close to the transition line where ferromagnetism 
disappears. Thus, they do not hint at the existence of a mixed ferromagnetic-glassy phase, 
as found in mean-field models,— in agreement with earlier T = numerical studies.— 1 ^ 

The FG transition line is slightly reentrant. Indeed, we find that p c = 0.7729(2) at T = 0.5 
and p c = 0.7705(2) at T = 1, which are definitely larger than p M = 0.76820(4), although 
they are quite close. This implies that there exists a small interval of the disorder parameter, 
around p ~ 0.77, showing three different phases when varying T: with increasing the tem- 
perature, the system goes from the low-temperature glassy phase with zero magnetization, 
to an intermediate ferromagnetic phase, and finally to the high-temperature paramagnetic 
phase. Correspondingly, it first undergoes a glassy-ferromagnetic transition with v = 0.96(2) 
and then a ferromagnetic-paramagnetic transition with v = 0.683(2). We mention that a 
slightly reentrant low-temperature transition line, where ferromagnetism disappears, also 
occurs in the phase diagram of the 2D ± J Ising model.— ^ 

The main features of the FG transition line are not expected to depend on the particular 
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FIG. 9: (Color online) Numerical results for the phase boundaries of the cubic-lattice ±J Ising 
model (HJ in the T-p plane. The dashed lines are interpolations of the data3. 

discrete bond distribution of the ±J Ising model, cf. Eq. ([2]). They should also apply to more 
general distributions with tunable disorder parameters, such as the Gaussian distribution 
reported in Eq. fl3]), and also to experimental spin glass systems with tunable disorder. 
We conclude showing Fig. M which reports all available numerical results for the phase 



.0 



boundaries of the cubic-lattice ± J Ising model ([I]) in the T-p plane, taken from Ref. |5| for 
the PF transition line, from Ref. |l9| for the multicritical point along the Nishimori (N) line 
T = 2/ln[p/(l — p)], from Ref. |6| for the data along the PG line, from this paper along the FG 
line, and from Ref. |l4| for the T = transition point. The dashed lines are interpolations 
of the data along the transition lines which satisfy the expected scaling behavior at the 
multi critica l point where they meet, controlled by the crossover exponent = 1.67(10), see 



Refs. 
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TABLE I: Results of combined fits of U4, U22, and R^ to Eq. ()Aip without scaling corrections, \ 2 
is the sum of the residuals in the fit and DOF is the number of degrees of freedom. Each column 
corresponds to results in which only data satisfying L > L min are included. R* = /r(0) is the 
value of the phenomenological coupling at the critical point. 



T = 0.5 

L min 4 6 8 10 

X 2 /B-OF 5594/289 567/229 203/169 79/109 

v 0.971(4) 0.964(5) 0.954(8) 1.00(2) 

p c 0.77230(1) 0.77275(2) 0.77284(3) 0.77281(5) 

#| 0.7453(2) 0.7564(3) 0.7592(6) 0.759(2) 

U% 1.3450(2) 1.3364(4) 1.3343(6) 1.334(2) 

Ufa 0.3046(2) 0.3045(3) 0.3057(6) 0.310(2) 

T = \ 

L min 4 6 8 10 

X 2 /DOF 10593/289 1842/229 365/169 98/109 

v 1.054(5) 0.995(5) 0.963(8) 0.982(20) 

p c 0.76819(2) 0.76920(2) 0.76975(3) 0.76994(5) 

Rl 0.6826(2) 0.7019(3) 0.7147(6) 0.7220(17) 

XJ% 1.3973(3) 1.3779(4) 1.3662(6) 1.3613(19) 

C/ 2 * 2 0.3094(3) 0.3020(4) 0.2977(5) 0.3013(17) 

Appendix A: Analysis details 

1. Phenomenological couplings 

In order to determine the exponent v and the critical parameter p c , we analyze the 
phenomenological couplings U4, U22, and R^ — £/L. In the critical limit each quantity R 
behaves as 

R(p, L) « f R [u p (p)L 1 ^] + u u (p)L- u g R [u p (p)L 1 /''], (Al) 

where the nonlinear scaling fields u p (p) and u w (p) are analytic functions of p. We have 
u p(Pc) = while, in general, we expect u u (p c ) 7^ 0. For both temperatures our data belong 
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to a small interval of values of p, so that we expect the approximations u p (p) ~ p — p c 
and Uoj(p) ~ u w (p c ) = a w to work well. To check it, we also performed fits assuming 
u p (p) = p — p c + k(p — p c ) 2 . We did not find any significant difference. 

We first analyze the results at T = 0.5. We perform combined fits of the three quantities 
to Eq. ( lAip without scaling corrections (we set a w = 0). If the scaling functions /r are 
approximated by fourth-order polynomials, we obtain the results reported in Table [H We 
report estimates for different L min : in each fit we only include data satisfying L > L min . 
Corrections are quite small and indeed the results corresponding to L min = 8 and L min = 10 
mostly agree within errors. We also perform fits that take into account scaling corrections. 
We fix uj, approximate gn{x) by a second-order polynomial, and repeat the fit for several 
values of uj between 1 and 5. If we perform a combined fit of U4 and R^ (we include all 
results with L > 4), the smallest x 2 /DOF (DOF is the number of degrees of freedom of 
the fit) is obtained for 3 < uj < 4 and one would estimate v = 0.96(1) and p c = 0.7729(1). 
If instead we use C/4, R%, and also U22 we obtain uj rj 2, v = 0.95(1), and p c = 0.7731(1). 
These results indicate that scaling corrections are quite small, and quite probably cannot 
be parametrized be a single correction term. Our best estimates of uj are simply effective 
exponents that parametrize the contributions of several different correction terms, which are 
all relevant for our small lattice sizes. 

If we compare all results, we end up with the estimates p c = 0.7729(2) and v = 0.96(2), 
reported in Eq. ([9]). For the phenomenological couplings at criticality, R* = /#(0), we obtain 
the estimates reported in Eq. JEJJ), i.e., R% = 0.764(6), t/J = 1.331(5) and £/ 2 * 2 = 0.305(2). 
The final estimates and their errors take into account the results of the fits with and without 
scaling corrections. 

The same analyses can be performed at T = 1. Combined fits to Eq. (1A1I) without scaling 
corrections give the results reported in Table fl] It is quite clear that scaling corrections at 
T = 1 are larger then those at T = 0.5. The goodness of the fit is worse and the fit results 
show systematic trends. It is however reassuring that they apparently converge towards the 
estimates ([9]) and ( flOl) . in agreement with universality. 

It is interesting to check whether scaling corrections can explain the differences which 
occur among the results for T = 1 reported in Tableland the results obtained at T = 0.5. 
Since the results for U^ at T = 1 are nonmonotonic as a function of L m i n , at least two 
correction terms must be included to explain the observed trend of the data. Therefore, 
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TABLE II: Estimates of the exponent 77 obtained by fits to Eq. (|A3p . where f x is approximated by 
a fourth-order polynomial and u(p) by a second-order polynomial. In each fit we only include the 
data which satisfy L > L m \ n . We fix v = 0.96(2) and the value of p c : p c = 0.7729(2) at T = 0.5 
and p c = 0.7705(2) at T = 1. 



T = 0.5 T = l 

AninrVDOF rj x 2 /DOF y 

4 516/94 -0.414(6) 62/94 -0.393(6) 

6 39/74 -0.400(8) 18/74 -0.389(9) 

8 22/54 -0.397(12) 16/54 -0.389(12) 

10 11/34 -0.398(16) 6/34 -0.390(16) 

the fit of the U22 data with a single scaling correction makes no sense. In any case, the 
estimate obtained for L min = 10 differs from the one reported in Eq. ffTUl) by one combined 
error bar, and therefore is in agreement with universality. We then perform combined fits 
of L/4 and £/L to Eq. ( 1A1I) . approximating Qr{x) by a second-order polynomial and fixing 
uj to several values between 0.5 and 1.5. The smallest x 2 /DOF is obtained for 0.8 < u < 
0.9. Correspondingly, we obtain p c = 0.7705(1), R£ = 0.765(10), and Ul = 1.32(1). The 
estimates of the phenomenological couplings at criticality are now in very good agreement 
with the estimates at T = 0.5. As for v we obtain v = 0.91(3), which is sligthly smaller 
than, but still consistent with the estimate at T = 0.5. If we fix v = 0.96(2) as obtained at 
T = 0.5, we find u = 0.95(10), p c = 0.7704(1), #| = 0.757(7), U* A = 1.326(6). 

These fits provide an estimate of p c at T = 1. We quote the estimate p c = 0.7705(2) 
already reported in Eq. f lT2|) . which satisfies the inequality p c > 0.7700, which one would 
obtain from the results reported in Table [H It is unclear how reliable our estimates of u are. 
In any case, they suggest a value close to 1. 

2. Magnetic susceptibility 

We analyze the magnetic susceptibility which should scale as 

X (P, L) = u h {pfL 2 ^f x [u p {p)L l/u l (A2) 
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TABLE III: Estimates of the exponent £• In each fit we only include the data which satisfy 
L > L min . We fix v = 0.96(2) and the value of p c : p c = 0.7729(2) at T = 0.5 and p c = 0.7705(2) 
at T = 1. 



T = 0.5 T = \ 

£minx7DOF C X 2 /DOF C 

4 88/98 -0.007(2) 161/98 -0.015(1) 

6 9/78 -0.003(2) 15/78 -0.009(2) 

8 8/58 -0.002(3) 3/58 -0.006(3) 

10 4/38 -0.005(5) 2/58 -0.005(5) 



where Uh(p) is related to the magnetic scaling field and is an analytic function of p; scaling 
corrections have been neglected. In order to determine r], we perform fits to 

In* = (2 - rj) InL + f x [(p - p c )L l ' v ] + u(p), (A3) 

where f x is approximated by a fourth-order polynomial and Uh(p) is normalized so that 
u(p = p c ) = 0. In this expression we have replaced u p (p) with p — p c . Inclusion of the 
second-order term does not change the quality of the fit and the results. Instead, even if 
the interval in p is small, the function Uh(p) cannot be approximated by a constant, hence 
u(p) cannot be set to zero. Indeed, the fits in which u{p) is approximated by a second- 
order polynomial have a x 2 /DOF which is significantly smaller than those in which we 
set u{p) = 0. For instance, for T = 0.5 and L min = 6 (we fix p c and u, see caption of 
Table HID, we have X 2 / D0F = 265/76 and 39/74 for the fit with u(p) = and the fit with 
a second-order polynomial, respectively. The results of the fits in which we fix v and p c are 
reported in Table [TT1 The results are very stable with L min and are completely consistent 
with universality. Note that, at variance with what is observed for the phenomenological 
couplings, corrections for T = 1 are apparently smaller than for T = 0.5. This may indicate 
the presence of several corrections which cancel out for our values of L. A conservative final 
estimate is r] — —0.39(2), already reported in Eq. (JTJ 
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3. Hyperscaling 

In order to study hyperscaling we consider the ratio 

H ee -^-. (A4) 

If hyperscaling holds, it should behave as 

H(p, L) = f h [u p (p)L l/v ] « f H [(p ~ Pc)L l '% (A5) 

where we have neglected scaling corrections. In order to allow for a possible hyperscaling 
violation we introduce a new exponent ( and assume that 

H(p,L)=L^f H [(p-p c )L 1 ^]. (A6) 

To determine ( we perform fits to 

\nH(p,L) = C\nL + f H [(p-p c )L 1 /^ (A7) 

where Jh{x) is approximated by a second-order polynomial. Fit results are reported in 
Table IIHI Here we fix v and p c to the values determined above. The quality of the fits is 
very good and scaling corrections are apparently small for both values of the temperature. 
The exponent ( is clearly compatible with zero, proving that hyperscaling is satisfied. More 
precisely, we obtain the bound |£| < 0.01, already reported in Eq. ffTS]) . 
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and 4> = 1.67. In Fig. [9] the FG line with T > 1 is approximated by the line of equation 28 
P ~ Pm = au 2(PiT)'t > , where a ~ 0.03 (which is fixed using our numerical estimate of p c at 
T = 1), U2<j>, T) = tanh/3 — 2p + 1, and pm is the position of the multicritical point. For T < 1 
we report the straight lines connecting the data at T = 1, 0.5 and T = 0.5, 0. Analogously, we 
proceed for the PG line, reporting the curve p — pu = cu2(p, T)^ for p > 0.7 and a straight line 
for 0.5 < p < 0.7. 
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